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Ti02 anatase doped with Co has been recently reported to exhibit room-temperature ferromag- 
netism. Ab initio study on substitutional Co doping, however, yielded much larger magnetic moment 
for Co than experiment. Our calculations based on density-functional theory show that the substi- 
tutional Co ions incorporated into Ti02 anatase tend to cluster and then the neighboring interstitial 
tetrahedral sites become energetically favorable for Co to reside, yielding a local environment more 
like C03O4 than CoTiOs. The interstitial Co destroys the spin-polarization of the surrounding 
substitutional Co but enhances the stability of the ferromagnetism significantly. In the absence 
of carriers, this room-temperature ferromagnetism can only be accounted for by superexchange 
interaction. 

PACS numbers: 75.50.Pp, 61.72.Ww, 75.30Hx 
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The discovery of ferromagnetism in IITV semiconduc- 
tors and the successful control of spin coherence of elec- 
trons injected from a magnetic semiconductor into a non- 
magnetic semiconductor suggests the possibility of har- 
nessing both charge and spin for new functionalities [? ? 
? ]. The limiting factor that represents a serious bottle- 
neck for their practical spintronic applications is the fact 
that both the observed ferromagnetism and the attrac- 
tive injection phenomena are essentially limited to low 
temperature. The development of new materials grows 
rapidly. Room-temperature ferromagnetism has been ob- 
served, for instance, in Mn doped ternary compounds 
such as CdGeP2[? ], ZnSnAs2[? ], and ZnGePsp ]. 

Using combinatorial pulsed-laser-deposition (PLD) 
molecular beam epitaxy (MBE) technique, Matsumoto 
et al.[l ] have fabricated Co-doped anatase thin films 
with different Co contents on SrTiO3(001) surface. At a 
concentration of 0.07, the film was reported to be ferro- 
magnetic at room-temperature with a magnetic moment 
of 0.32 per Co. In a more recent work. Chambers 
et aL[? ] employed oxygen-plasma-assisted (OPA) MBE 
to grow Co-doped anatase films on the same substrate. 
They observed a significantly larger magnetic moment for 
Co ions, i.e., 1.26 ^b, which is much closer to the value 
of the low spin state of Co. From the comparison of 
the surface-sensitive X-ray photoemission spectrascopy 
(XPS) and bulk-sensitive X-ray absorption and emis- 
sion spectroscopy (XAS) of Co-doped Ti02 anatase and 
various other reference compounds such as Co, C02O3, 
C03O4, and CoTiOa, Co ions in anatase were judged to 
be substitutional and in the -1-2 formal oxidation state. 
Also, they found the Co distribution, and therefore the 
magnetic properties depend critically on the growth con- 
dition. 

To gain a microscopic understanding of the magnetism 
of this noval materials. Park et aL[? ] carried out 
a systematic computational study on Tio. 9325-^0.062502 
(M=Co, Mn, Fe, and Ni) using the linearized muffin-tin 



orbital (LMTO) method both with local-spin-density ap- 
proximation (LSDA) and with LSDA -I- U incorporating 
the on-site Coulomb correlation interaction U . Their cal- 
culations showed that the spin magnetic moment of Co 
was about l/xs, indicating a low spin state. But they also 
obtained an orbital moment for Co as large as 0.9/iB, so 
the total magnetic moment per Co was 1.90/iB, about 
six times as high as the PLD-MBE result [? ], or, 50% 
higher than the OPA-MBE experiment [? ]. This dis- 
crepancy is obviously too large to be acceptable. Recent 
density-functional computations on various diluted mag- 
netic semiconductors such as Mn-doped GaAs[? ? ] and 
ZnGeP2[? ], for instance, have reproduced a magnetic 
moment of Mn in good agreement with experiment [? ? 
] . Although LSDA is not appropriate to describe strongly 
correlated systems, LSDA-I-C/ has proved quite success- 
ful to yield a valid magnetic structure for such systems [? 
]. Thus, the marked disagreement between experiment 
and theory casts some doubt on the understanding of 
the magnetism and its connection to the structure of Co- 
doped anatase. Motivated by the recent discovery that 
defects could play a key role in the magnetic structure of 
the diluted magnetic semiconductors [? ? ], we conjec- 
ture that clustering of substitutional Co and/or intersti- 
tial Co, if exists, might be responsible for the reduction 
of Co magnetization. 

In an attempt to clarify this point, we have performed 
ab initio density-functional calculations to study the 
structural, electronic, and magnetic properties of Co- 
doped Ti02 anatase. We show that the substitutional 
Co ions tend to cluster, but this has no remarkable effect 
on the magnetization of Co. We also find that a (seri- 
ously flattened) tetrahedral interstitial site will become 
energetically favorable for a Co atom (in reference to bulk 
Co) when two of its neighboring substitutional (octahe- 
dral) sites are occupied by Co. The interstitial Co is 
also in the low spin state (1.04^^), but it exerts strong 
detrimental effect on the magnetization of its neighbor- 
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ing Co —0.04/1^) and therefore reduces the 

magnetic moment Co in this magnetic semiconductor. 
While substitutional Co ions turn anatase half-metallic, 
the interstitial ones bring the material back to a semicon- 
ductor. In absence of carriers, the exchange interaction 
stablizing the ferrimagnetism is of superexchange type. 
In the following, we give a brief outline of our theoretical 
procedure, a discussion of the results based on bonding 
characteristics and its relevance to oin- understanding of 
the magnetism in this Co-doped oxide. 

The calculations were done using a spin-polarized ver- 
sion of the Vienna ab initio Simulation Package (VASP)[? 
]. Fully nonlocal Vanderbilt-type ultrasoft pseudopoten- 
tials[? ] were used to represent the electron-ion inter- 
action. Exchange correlation interactions were described 
by the Perdew-Wang 1991 generalized gradient approx- 
imation (GGA)[? ]. Wave functions and augmentation 
charges were expanded in terms of plane waves with an 
energy cutoff of 396 and 700 eV, respectively. The Bril- 
louin zone integration was performed within Monkhorst- 
Pack scheme [? ] using a (2x2x2) k mesh for geometry 
optimization and a (4x4x4) mesh for plotting the den- 
sity of states(DOS). For each system, the geometry was 
relaxed until the atomic forces were smaller than 0.03 
eV/A. 

To gain a confidence for applying pscudopotentials to 
such an oxide with open structure, we first conducted 
a comparitive study of the structural properties of non- 
doped anatase (see Figure la) using both VASP and the 
highly precise all-electron full-potential linearized aug- 
mented plane wave (FLAPW) method[? ] with WIEN97 
implementation[? ]. FLAPW is considered to be more 
accurate but also more computationally expensive. The 
same GGA functional and k mesh (8x8x4) were used 
in both calculations. The optimized structural param- 
eters for anatase are listed in Table I. Also listed are 
experimental values [? ] and a previous FLAPW -I- LDA 
result [? ]. It is evident that VASP agrees fairly well with 
WIEN97, and that GGA results are more consistent with 
experiment than LDA. 

For Co-doped anatase, we mainly dealt with a supercell 
containing 16 primitive unit cells (2a x 2o x c, sec Figure 
lb) with one or two Ti atoms (site Bl, or, sites Bl and 
B2) replaced by Co and with/without an interstitial Co 
(site A). In view of the fact that the calculated c/a value, 
2.46, of a hyperthetical C0O2 in anatase structure is close 
to that of Ti02 anatase, we fixed c/a at the latter value 
when optimizing the volume of Co-doped anatase to save 
computational effort. To see whether substitutional Co 
(Co'^) cluster or not, we first minimized the free energy 
of the cell Tii5Co;Q^032. Upon Co-doping, the volume 
of such a unit cell contracted by 1.1%, in accordance 
with the fact that Co has a smaller covalent radius than 
Ti[? ]. Then we doubled this cell along x direction and 
brought the two Co'^ together to form Ti3oCo§;^Co§2064 
(4o X 2a X c). We assumed that there was no volume 
change upon clustering, and optimized only the internal 
freedoms for the 96 atom cell. We find the latter config- 



uration is more stable by 0.14eV/Co, indicating that Co 
clustering does occur. 

The calculated formation energy of an interstitial Co 

(Co^) as a neutral defect, £;(Tii6Co^032)--B(Tii6032)- 
£^(Co), is -1-2.07 eV, implying a sole Co^ is highly unsta- 
ble in reference to bulk Co. This can be understood in the 
bond-order band-strength picture. Ti:0 + Co:0 is much 
weaker than 0:Ti:0, as is evident from the comparison of 
the formation heat of Ti02, TiO, and CoO[? ]. Nonethe- 
less, with the presence of Co^j^, the formation energy of 
Co^ i.e., £;(Tii5Co|iCo;^032)-^(Tii5Co|i032)-S(Co), 
drops drastically to -f0.18 eV. According to our calcu- 
lations, the formation heat of C0O2 in a hypothetical 
anatase structure is —3.66 eV, slightly higher than two 
times of that of CoO (—1.59 eV). This means that Co:0 
-|- Co:0 is only slightly weaker than 0:Co:0 , thus Co^ 
becomes much less unstable. Further calculations show 
that when both Bl and B2 site arc occupied by Co upon 
clustering, site A becomes favorable for a Co^ with a 
formation energy of —0.49 eV. Once again, this can be 
explained by the higher formation heat of C03O4 (—9.57 
eV) than two times of that of C0O2 (—3.66 eV). Our 
calculations thus reveal a 2Co'^-|-Co''^ local structure in 
Co-doped Ti02 anatase. 

Table II displays the calculated Co-0 bond length in 
different local environments. The overall feature is that 
the Co-0 bond is shorter than Ti-0 bond. As mentioned 
above, this is because both octahedral and tetrahedral 
Co ions have a radius smaller than octahedral Ti ions. 
Co-0 bond lengths along the y+ and y— directions, dy+ 
and dy- , do not change upon clustering in the case that 
Co(Bl) and Co(B2) are in the same x — z plane; whereas 
along the x and z directions the bond lengths are con- 
tracted by 0.02 - 0.05 A. The presence of an interstitial 
Co pushes the surrounding atoms a bit further apart and 
increases the Co-0 bond length in the case of Co(Bl) 
and Co(B2). The local lattice distortion induced by Co"^ 
or Co^, however, is rather small. It is worth noting that 
the volume effect on chemical bonding is only marginal. 
A sole Co''^ yields a volume expansion of about 1.0% to 
Tii6032, and the formation heat of Tii6Co^032 at the 
volume of Tii6032 is only 0.05 eV lower than the opti- 
mized value. 

To examine the effect of clustering on magnetism, 
we performed calculations on both intra-cell FM 
and intra-cell AFM alignments of the Co pair in a 
Tii4Cog]^Co32032 unit cell (For simplicity, we denote 
this cell as Tii4Co2032 in the following) while keeping 
the inter-cell coupling ferromagnetic. The FM phase is 
found to be lower in total energy by 0.09 eV than the 
AFM phase. When the two Co ions are far apart, this 
energy difference is 0.08 eV[? ]. Apparently, clustering 
has negligible effect on the spin-polarization of this mag- 
netic oxide. In the case of 2Co"''+Co^, we checked all the 
three possible intra-cell spin alignments of the three Co 
ions. All of them converge to the same final state, with 
both Co^ and Co"^ in the low spin state. We then double 
the unit cell along x direction and computed the inter- 
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cell AFM coupling between two 2Co"^+Co^ complexes. 
It is found to be 0.35 eV higher in energy than the FM 
coupling. Obviously, with a much enhanced FM-AFM 
energy difference, the formation of 200*^+00^ will raise 
remarkably the Curie temperature of Co-doped anatase, 
in accordancic with the observed room-temperature fer- 
romagnetism[? ? ]. 

The calculated total spin magnetic moment in a unit 
cell of Tii5Coi032, Tii4Co2032, and Tii4Co3032, and 
that in the atomic sphere of Co are listed in Table III. 
It is seen that in all three cases, Co ions are in the low 
spin state, regardless to their oxidation states. The spin 
moment of Co"'' does not change much upon clustering. 
Co^, on the other hand, behaves rather differently. When 
coupled with Co'^, Co^ kills almost entirely the spin mo- 
ments of Co'^. As a result, the average spin magnetic mo- 
ment of Co is reduced to 0.32 ^b, less than a half of the 
non-Co^ case. This means that if there are a remarkable 
amount of 2Co'^-|-Co^ complexes in Co-doped anatase, 
the magnetic moment per Co atom would be reduced 
significantly and this could explain the disagreement be- 
tween experiments [? ? ] and the previous first-principles 
calculations [? ]. We note that the oxygen-rich condition 
during OPA-MBE growth in Chambers et alJs work[? 
] prevented the formation of extensive 2Co'^-|-Co^ com- 
plexes, hence a much higher magnetic moment for Co 
than that in the film grown by Matsumoto et aL[? ] 
with PLD-MBE; first-principles calculations [? ] gave a 
even higher magnetic moment for Co'^ , as a result of the 
absence of Co^. In oxygen-poor condition, on the other 
hand, the oxygen vacancies help the diffusion of the Co 
ions, resulting in the formation of Co metal clusters [? ]. 

Figure 2 displays the calculated density of states 
for non-doped and various cases of Co-doped anatase. 
Clearly, the Co states are mainly located in the energy- 
gap region and the O and Ti states are not much af- 
fected by Co doping. A comparison of Tii5Coi032 and 
Tii4Co2032 indicates clustering does not change the half- 
metallic electronic stnicture, and as a consequence the 
magnetic interaction between Co ions remains to be 
carrier-mediated Zener-like. [? ] Sole Co^ ions turn the 
material into a metallic system with a highly lifted fermi 
level. Very interestingly, our calculations show a complex 
of 2Co^-|-Co^ brings the system back to a semiconductor, 
in consistence with the experimental observation [? ? ]. 

To gain a deeper insight into the electronic structure 
of Co, we compare d DOS of Co in anatase with those 
of Co in other Co oxides with different oxidation states 
such as C0O2 (+4), CoTiOg {+2), and C03O4 (+3 and 
+2). C0O2, the end member of LiCo02, has recently 
been isolated [? ]. It seems to be of CdCl2 type with 
a monoclinic distortion, but positions of O remain to be 
determined. In this regard, we believe a study on a hypo- 
thetical C0O2 in anatase structure is meaningful in illus- 
trating the electronic structure of Co*+. As mentioned 
above, we have optimized all of its structural parameters. 
For CoTiOs and C03O4, experimental lattice constants[? 
] were adopted and only internal freedoms were opti- 

FIG. 1: Calculational unit cells of non-doped (panel a) and 
Co-doped Ti02 anatase (panel b). Large and small circles 
represent O and Ti (Co), respectively. 



FIG. 2: Total density of states for non-doped and various 
cases of Co-doped Ti02 anatase. Dotted vertical lines denote 
the Fermi energy. 



mizcd. The calculated d DOS of Co in C0O2, CoTi03, 
and C03O4, and that in Co-doped anatase are plotted in 
Figure 3. We observe that although much more localized, 
d DOS of Co^ in Tii5Coi032 and Tii4Co2 032 resembles 
that of Co in C0O2 in both exchange splitting and crystal 
field splitting, suggesting they are in +4 oxidation state, 
as was also proposed in Ref.[9]. In an octahedral field, t2g 
is lower than eg and thus Co^+ will have a spin magnetic 
moment of I/xb- On the other hand, Co''^ in Tii6Coi032 
reminds us of Co in CoTiOs. This resemblance reveals 
Co^ is in +2 oxidation state in a Tii6032 environment. 
In a normal tctrahcdral field, Cg level is lower than t2g. 
Co^ , nevertheless, is in a seriously flattened tetrahedral 
field, thus the crystal field splitting states is comparable 
to the exchange splitting, yielding an intermediate spin 
state for Co'^. As for the combination 2Co"^+Co^, an 
even better similarity can be found between Co'^ (Co''^) 
in anatase and Co^ (Co^) in C03O4. It thus provides a 
strong evidence that Co'^ and Co^ in such a combination 
are in -1-3 and +2 oxidation states, respectively. Simi- 
lar to Co^ in C03O4, Co'^ in an octahedral crystal field 
has a vanishing spin-polarization, and the contribution 
to magnetism comes only from Co^. 

To conclude, our prediction of the formation of 
2Co'^-|-Co-^ resolves the disagreement between experi- 
ment and theory on the magnetic moment of Co in Co- 
doped TiO^ anatase. The insulating property is repro- 
duced without employing an on-site Coulomb correla- 
tion interaction U . Supcrexchange, rather than carrier- 
mediated Zener-like magnetic interaction, should be re- 
sponsible for room-temperature ferromagnctism. The 
novel strucure discoverd in this work may also exist 
in other transition elements such as Mn- and Fe-doped 
anatase and may thus resolve the discrepancy between 
experiment [? ] and theory [? ]. Finally, we want to 
note that the formation of more complicated (Co*^, Co^) 
combinations is also possible depending on the growth 
condition, but the formation probability should be lower 
than 2Co'^-FCo'f. 
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FIG. 3: The calculated density of d states of Co in C0O2, 
CoTiOs, C03O4, and various cases of Co-doped Ti02 anatase. 
Dotted vertical lines denote the Fermi energy. 



TABLE I: Optimized structural parameters for anatase Ti02 : 
A comparison of all-electron full-potential method (WIEN97), 
pseudopotential method (VASP), experiment [? ], and a pre- 
vious FLAPW LDA study [? ]. 





WIEN97 
(GGA) 


VASP 
(GGA) 


Expt. 


FLAPW 
(LDA) 


a (A) 


3.826 


3.825 


3.782 


3.692 


c(A) 


9.706 


9.678 


9.502 


9.471 


c/a 


2.537 


2.530 


2.512 


2.566 


u 


0.207 


0.207 


0.207 


0.206 



TABLE III: Claculated spin magnetic moment (/is) in the 

atomic sphere of Co (radius = 1.10 A) and that for a whole 
unit cell. 

Unit cell Co(Bl) Co(B2) Co (A) Total 

Tii5Coi032 0.73 - - 1.00 

Tii4Co2 0:i2 0.69 0.69 - 2.00 

Tii4C():!():!2 "O.Ol -0.1)7 i.Ol i.OO 



TABLE II: Optimized Co-0 bond lengths (A) in Co doped 
anatase Ti02. It is denoted as dx+ {dy+, dz+) when O is in 
the x+ direction of Co, and as dx- {dy-, dz-) when O is in 
the X- direction of Co. Co(Bl), Co(B2), and Co(A) denote 



Co ions 


occupying 


sit(-s Bl 


. B2. and 


A (sec Fi.<j 


,urc lb). 


dco-o 


Tii5Coi032 Tii4 


C02O32 


Tii4 


C03O32 


Co(Bl) 


Co(Bl) Co(B2) Co(Bl) Co(B2) Co(A) 


dx+ 


1.91 


1.85 


1.93 


1.93 


1.98 1.92 


dx- 


1.91 


1.93 


1.85 


1.88 


1.93 1.91 


dy+ 


1.91 


1.91 


1.91 


1.96 


1.97 1.91 


dy- 


1.91 


1.91 


1.91 


1.90 


1.93 1.85 


dz+ 


1.89 


1.84 


1.95 


1.93 


1.98 


dz- 


1.89 


1.94 


1.84 


1.96 


1.98 




-4 -2 2 4 6 8 

Energy (eV) 



